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AXIAL COMPRESSOR DESIGN FOR 
AIRCRAFT ENGINES. 


By Norman C. BrEppy. 


THE axial flow compressor is based on the conversion of kinetic 
energy to pressure energy through the medium of rows of moving 
blades and stator blades which are designed on an aerodynamic 
basis. As such, the axial compressor provides a smoother gas flow 
than the centrifugal type and greater pressure ratios are obtainable. 
From the standpoint of energy transfer the compressor is simply 
a reversed turbine but there the similarity ends. 

The principle of the axial flow compressor is by no means a 
new invention. C. A. Parsons filed a patent in 1884 which described 
the use of his reaction turbine as a compressor ; however, it is only 
during recent years that the difference between an accelerating 
flow and a decelerating flow has been fully understood. 


In all flow phenomena there exists an associated boundary 
layer in contact with the containing walls, which, because of 
friction, flows at a decreased velocity. There is a natural tendency 
for the gas to break away from the walls in a diffusing passage 
and cause turbulence, or in extreme circumstances result in a 
backflow. | Unless the decelerating process is carefully controlled 
by aerodynamic theories the boundary layer can become so thick, 
that the main flow is contained within a channel of constant cross 
section illustrated in Fig. 1. In such a case P, would be only 


P, 22 a A LS aS AaB | 


7999 9,9 
2 oo = 


Fig. 1—Inefficient Diffusion Process. 
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slightly greater than P, and therefore practically no velocity change 
has occurred. During recent years a tremendous amount of data 
has been obtained by the study of flow past a series of aerofoils 
known as a cascade, which represents the flow problem of a com- 
pressor with a greater degree of accuracy than by using a single 


aerofoil. 

Nomenclature. 

U = Tangential velocity of blades in feet per sec. 

Va = Axial rey of air through the compressor in feet per 
second. 

V, = Relative velocity of air in feet per sec. 

oc = Air angle. 

Vw = Whirl cémponents of the absolute air velocities in ft. {sec. 

Ww = Work-done in ft. Ibs. ’ 

Q = Work done factor. 

ATs = Stage total head temperature rise. 

Rst = Stage total head pressure ratio. 

Ti: = Inlet total head temperature. 

Nst = Total head isentropic efficiency. 

AT. = Static temperature rise in rotor. 

AT, = Static temperature rise in stator. 

A = Degree of reaction. 

T = Static temperature. 

“A = Area in sq. ft. 

@ = Density of air in Jb./cu. ft. 

L = Lift in lbs. 

D = Drag in Ibs. 

C.p = Lift co-efficient. 

C,p = Drag co-efficient. 

Pit = Inlet total head pressure 

Pot = Outlet total head pressure 

P; = Inlet static pressure Ibs. per sq. in. 

P, = Outlet static. pressure 

AP = _ Static pressure rise 

F = Axial force in lbs. 

Q = Mass flow in Ib. per sec. 

y = Ratio of specific heats 7 

T2,1 = Outlet total head temp. for isentropic compression. 


General Aerofoil Theory. 


Fig. 2 (a) represents an aerofoil supported in a wind tunnel 
in a fixed position with reference to the air stream; 1.¢., the flow 
is parallel to the aerofoil. In this case the air is divided around 


COMPRESSOR DESIGN FOR AIRCRAFT ENGINES 5 


the section, separating at the leading edge and joining again at 
the trailing edge. With a well-designed aerofoil there exists no 
turbulence and the flow is streamlined. 

Now let the aerofoil be set at an angle 7 degrees, which is termed 
the angle of incidence, to the air flow, as shown at Fig. 2 (b). In 
this case a greater disturbance is set up and the flow will be changed 
although at a point well behind the trailing edge the airstream will 
resume its normal flow. These local deflections can only occur 
when the blade exerts a force on the air, thus by Newton’s law an 
equal and opposite reaction takes place on the blade through the 


Fig. 2—Air Flow Past an Aerofoil. 
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medium of the air. Such forces can only appear in the shape of 
pressure exerted by the air stream on the blade. 


By the Bernoulli equation it follows that with local pressure 
changes along the blade profile there exists velocity changes. 
Examination of the streamlines over the top of the aerofoil shows 
that they tend to converge upon each other, thus indicating an 
increase in velocity, with a corresponding drop in pressure. On 
the lower side of the aerofoil the streamlines separate indicating 
a reduced velocity with an increase in pressure. The pressure 
distribution around such an aerofoil is shown in Fig. 2 (c). 


The vectorial sum of these pressures will result in a resultant 
force which can be resolved into L, a lift component and D, a drag 
component. It is possible to measure the lift and drag force for 
= Aes of aerofail sections, airflow and incidence angle, etc., 

1g. 9. : . : 


Now lift L = ae ee (1) 
2g 
1 
Drag D = ro O..V*.CScA (2) 


From (1) and (2) has been evolved two co-efficients C, and C, 
relating velocity, area, density and lift or drag forces. Curves 
can be plotted from wind tunnel tests which can then be employed 


- in calculations relating to the particular aerofoil. The aerodynamics 


of flow throughout a compressor can be analysed in two ways; 
(1) as an isolated aerofoil, and (2) as a cascade of aerofoils. The 
second method is now the accepted one and Fig. 4 illustrates the 
cascade nomenclature for compressor blading. 


Nomenclature, 
8, = Blade inlet angle. 
8B, = Blade outlet angle. 
6 = Blade camber angle = B, ~ Bz. 
¢ = Stagger angle. 
S = Pitch. 
C = Chord. 
a = Distance of point of max. camber from the leading 
edge. 
o, = Air inlet angle. 
oc, = Air outlet angle. 
V, = Air inlet velocity. 
V, = Air outlet velocity. 
t = Incidence angle = oc, -§8,: 
5 = Deviation angle = oc,—8,. 
e = Deflection = ,— oy. 


COMPRESSOR DESIGN FOR AIRCRAFT ENGINES 7 
Lift and Drag Co-efficients of a Cascade of Blades. 


In practice it is usual to relate the lift and drag of the cascade 
to the vector mean (Vm) of the inlet and outlet gas velocities. The 
axial velocity is assumed to be constant throughout the stage and 
under such circumstances the velocity diagram of the cascade is 
shown in Fig. 5. ; 

The static pressure rise across the blades is given by :— 

AP = P,-P, = (P-$e Vs!) = (Pu-de Vs 
AP =4e@ (V,2-V2) —@ = 4 Va? (tan®oc, —tan®ec,) -@ (3) 

Due to the fact that the density change is negligible the in- 

compressible flow formula is adopted. 


RESULTANT FORCE 


L (ut) 


ANGLE OF INCIDENCE 
4 


b(DRAC) AIR FLOW 


LIFT (LBS) 
DRAG (LAS) 


-t0 0 0 IS 20 


ANGLE OF INCIDENCE (7) 
Fig. 3—Forces Acting on an Aerofoil. 
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The axial force F = S.g.V,x change in velocity component 
along the cascade = S.g.V,? (tana, —tanc,) (4 
Now Vm = Va sec Om (5) 
and tan oc, = } (tanoc,+tano,) (6) 
D and L being the drag and lift forces along and perpendicular 
to the vector mean velocity 
D = 4e@Vn?.C.C, = Fsin em-S. AP. cos xm (7) 
From equations (3) and (4) we get :— : 
4@.Vm?.C.C,p = S.Q. Va? (tan o,—tan o,) sino, —S.4e.V." 
(tan?oc, —tan%oc,) cos Om+S.COS Sm.w 
and since (tan*oc, ~tan?oc,) = (tan o¢,—tan oc.) (tanoc, +tanoy) 
= 2 (tan oc, -tan o,) tan Cp 


S w cos Gm S w  cos*op 
Cop = FoR ga ea I ae ee eee 
C te ° Vm? C te Va 
S- w '  costacy 


= _——_ .—— (8) 
Cc s0V? cos*oc, 
By resolving forces perpendicular to the vector mean Vian 
L=4e@.Vn?.C.C, = Fcoos om+S. APsin <q _ 
Therefore $@ Vn?.C.C, = S.@. Va? (tan oc, —tan o,) cos cm —w 
S.sin Cp 


S 
“ G=2. Cc (tan oc,—tan OC.) CoS Sm—C,p. tan Ca (9) 

The term C,p tan oc, in equation (9) is so small that it can be 
neglected hence :— 


Cc. = 2. = (tan oc, —tan oC.) cos Om (10) 


Gascade Performance. 


Tests are now made ona cascade of blades of known shape at 
a given stagger angle with a range of incidence and flow velocities. 
Measurements are taken of the inlet and outlet velocities (V, & V,) 
along with the air angles (oc, and oc,) and the pressure difference 
across the blade row (P,—P,). From this data the following 
quantities can be calculated :— 
Cm (Equation (6) ). 
Vim (Equation (5) ). 


€ a Oy = Co. 
t+ = «,-f,. 
3) 


= ,-f,. 

AP (Equation (3) ) and pressure loss co-efficient AP/$,.g.V,”. 
Cyp (Equation (8) ). 

C, (Equation (10) ). 
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Fig. 5—-Velocity Diagram (Constant Axial Velocity). 
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Fig. 6 illustrates the performance of a compressor cascade. 
It will be seen that a value of deflection can be selected which is 
the most suitable for the cascade being considered. The ultimate 
object of the cascade is to turn the air through the largest angle 
possible, with a minimum loss, hence the selected deflection will 
be near the maximum point of the curve. 


C 
‘In order to obtain a high efficiency the ratio + should 
DP 

be at a maximum and in this particular example it is about 50. 
At this figure the deflection is 80% of the maximum deflection. 
Successive tests carried out on other cascades show similar character- 
istics ; thus it has beceme a standard practice to adopt a “design 
point’”’ on a cascade where the deflection is 80% of the maximum 
obtainable. Design’ conditions are usually denoted by a star, 
t.e., €* = 0:8 (oc, -oc,) max. (11) 

For blades of the same pitch/chord ratio and camber, equal 
values of «* are obtained when the corresponding outlet angles 
oc, are also equal. Although « * at any given oc, decreases as 
the pitch S increases, the variation with the camber angle @ is only 
very small and over the range 20-40 degrees, can be disregarded. 
Therefore design deflection ¢ * is a function of S/C and a, See 
Fig. 6. 

The air outlet angle oc, for a given blade outlet angle f, is 
measured by the deviation 6 which can be estimated from a number 
of empirical equations. Constants Rule is widely used :— 


) 6g 3 ; 12 
= mo (c) | (72) 
For blades of circular or parabolic camber 
2a,” ee 
- 023 (=) +01 (=) 13 
m c). + 50) (13) 
When a circular arc camber line is used (13) becomes :— 
oc 
m = 023401 ( 2) 14 
50 (14) 


Work Done Por Stage. 


Fig. 7 represents one stage in an axial flow compressor. Air 
approaches the rotor blades with an absolute velocity V, and at 
an angle cc, in combination with a peripheral velocity U. The 
relative velocity will be V, at an angle oc, as shown by the velocity 
triangle. The air then passes through the diverging passages 
which increases its absolute velocity to V, at an angle oc, Due 
to diffusion V, < V, and thus a pressure rise would have occurred 
in the rotor. The peripheral velocity U in combination with V,, 


DEFLECTION (DEGREES 


LIFT COEFFICIENT 


DRAG COEFFICIENT 
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IRCIDENCE 1 DEGREES 


Fig. 6—Performance of a Compressor Cascade. 
For @ = 38°, A/C = 0-4, S/C = 0-94, B, = 14° (from A. R. Howell). 
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gives the absolute velocity V, and outlet angle oc, at exit from the 
rotor. 


From the geometry of the velocity triangles two basic equations 
are obtained :— 


U 

vy, = tanoc, + tan oy (15) 

U 

yo tan cc, + tan o, (16) 
a 


The work done on the air in the rotor blades per lb. of air per 
second is equal to the change in whirl momentum times the distance 
moved per second. . 


WU 28 U 
W = Boe een eo) = Pia 
(tan oc, —tan o,) (17) 
The energy imparted to the gas will increase the total head 
temperature a theoretical amount. 


UV 
AT, = — > (tan oc, —tan oc (18) 
g : J . Cs ( 1 2) 
For J = 1400 ft. lb. per C_H.U. and C, = 0-24 equation (18) reduces 
5 to 
U Va 
AT, = 0-923 Greg. ( 10 ) (tan oc, -—tan oc,) (19) 


The actual stage temperature rise will be somewhat less than 
this, due to the fact that the axial velocity V is not constant along 
the blades. The increase in total head temperature is only about 
85% to 90% of equation (19) thus :-— 

U Va 
ATs: = Q.0-923 ered, Ga (tan oc,—tan oc,) (20) 
where {2 = work done factor. 
st. ATst yl 
The stage total-head pressure ratio Rst = (1 a ae 
1 


(21) 


where T,¢ inlet total head temperature. 
mst = total head isentropic efficiency. 


Degree of Reaction. 


Reaction is the ratio static temperature rise in the rotor to that 
in the whole stage. Under isentropic conditions this is equal to 
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Fig. 7—Axial Compressor Stage. 
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the corresponding ratio of static pressure rises. The choice of 
the degree of reaction in compressor design is an important one 
and an expression can be evolved for a certain degree of reaction 
in terms of velocities and air angles. The simplest case is one 
where it is assumed that V4 = V,, 7.¢., the gas leaves the stage 
with the same velocity at which it enters it, hence ATs = ATst. 


Putting AT, = temperature rise in rotor. 

; AT, = temperature rise in stator. 

Then W = Cy.J (AT,+AT,) = Cp.J. AT; 
U.Va 


(tan oc, —tan 4) (22) 


e 


U.Na 
= (tan ocs—tan oc,). 


Vo = Va sec x, 
*, Vo?-V,? = Va? (tan?oc, — tan?oc,) (23) 
Vs = Va sec a, 


a 
By substitution T,-T, = ——— (tan oc,—-tan oc, 
> a—4y g.J.c% ( 3 ) 
| mle tan oc,—tan oc,) 24) 
— 2U ( 3 Lo ( J 
AT, AT, 
The degree of reaction, 7.¢., = 
- ‘ AT; AT, + AT, 


_ U.Va (tan ,—tan o,) — $Va?* (tan*oc, — tan®oc,) 
U.Va (tan oc, —tan oc,) 
Va 


U 


tan oc, + tan o, 
tan oc, + tan o, 
In the case of 50% reaction blading where the inlet and exit 


air angles of the stator are equal to the inlet and exit angles of the 
rotor. 


oc, = o,and oa, = &, .. from (25) A = 1-$ 


Vortex Flow. 


Consideration has so far only been given to the velocity triangles 
at one particular radius. The peripheral speed U varies proportion- 
ately to the radius and the variation of the other velocity vectors 
will depend on the design. The velocity of air at any point during 


oe 


(tan oc, + tan oc,) 


= 1-4 (25) 
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the flow through the compressor will possess a component in the 
peripheral as well as the axial direction, therefore the flow will be 
a spiral in the direction of rotation, otherwise known as vortex 
in form. Vortex theory is based on the fact that radial pressure 
forces act on the air elements in order to obtain the necessary 
centripetal acceleration in association with the whirl velocity. 
That is, the whirl velocity varies inversely as the radius which in 
equation form :— 


constant 
Vo = oS 


7 


Stage Performance and Pressure Rise. 


In addition to the drag co-efficient already mentioned, two 
additional factors must be taken into account. They are secondary 
losses, attributed to additional drag in the compressor annulus, 
tip clearance and trailing vortices. Tip clearance has a big 
influence on such losses, thus should be kept as small as possible, 
usually in the region of 1% to 2% of the total blade height. The 
additional drag co-efficient arising from secondary losses has been 
approximated by the following empirical formula :— 

Ci. = OOISC? .. (26) 
C, being the lift Co-efficient given in equation (10). 
_ The annulus drag loss has also been quoted as an empirical 
formula :— : 
C,. = 0:020 s/h (27) 
where s = pitch. 
h- = height of blades. 

C, = Cop en Cys + Con 

Equation (8) was derived for the straight cascade and will apply 
equally as well for the annular case. Substituting C, for Cp 

S w 3 Om 
Cu ie: 5, w _ Losec (28) 
C 4$@.V,? cose, 


The overall pressure loss for the stage —_ 
3@.Vi 
Cy 
5 S_ _cos*&m , (29) 
Cc cos?0c, 
If the compression in the blade row takes place according to the 


P 
polytropic law 2" = C 
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the efficiency of the process is 


Tit =Tet 
Tot - Tye 


‘Ca ae 
Pad 
——_—_—_—_— (30) 


Putting the loss w equal to zero and substituting in equation (3) 
the theoretical pressure rise 


AP! = 4@.V,? (tan? oc, —tan? o.) 


e 


API. : 


4@. Va? sec? a, 


ASP? 
ze@.Vi 


The efficiency of the 


4@.V,? (sec? oc, - 


sec? 0») (31) 
5 pe. sec? oy 
sec? oc, 
cos? oc, 
oe cos* oy, 
i= 
blade row 7, = a = 1|- 
Gt 2: 
$@.Vi? 
—~ APL (32) 
z@.V, 


With 50% reaction blading this formula will hold good for both 
stator and rotor blade rows. 


Also 7, 


P,—-P, 
P,t-P, 


The stage efficiency 7, is the ratio isentropic temperature rise/ 
actual temperature rise. 


P. 
P, 


ATs 


P,} 
P, 


= [1+ 


y-1 


E i nsATs ic 
oF. 


being the temperature rise in the rotor for 50% reaction. 


Po 

— _-} 
ATs Vy. P, 
oT, Fe 1 Pt 
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nsAtTs na 
14 SS— la - 
L * “er, I? 


ATs “Ly. 
1 = —i1 
L * “OT, ia 


” 1 ATs 
= 1- ——~.—(l- 34 
n[ y-l 47, ( 1) | 64) 
In present day power plants AT, would be approximately 35°C. 
and T, about 288°K. at the first stage thus the second term in the 
bracket can be disregarded giving :— 


Qe = 7s 
In the case of the half vortex case an approximate stage efficiency 
is evolved by taking the arithmetic mean of the efficiencies of the 
two blade rows. 
ns = 4 (yn, Rotor + 7s Stator) (35) 
When the degree of reaction is considerably removed from 50% 
reaction a more accurate expression for the stage efficiency would 
be :— 
Ns = A-] Rotor + (1— A). 7s Stator (36) 
The degree of reaction a being obtained from equation (25). 


Characteristics at High Speed. 


The critical condition at any point in a blade passage is when 
the velocity of the air exceeds that of the local speed of sound, 
which causes shock waves and a loss in total head pressure. The 
velocity of sound in air increases with rising temperatures thus 
the Mach numbers in a compressor will decrease progressively 
through the stages. On this basis the early stages are the most 
critical and the worst condition is when the compressor is operating 
at altitude. . 


Va 


(y.g-R.T)) 
The performance falls after the Mach number reaches 0-7 and 
gradually increases until the flow chokes at the minimum cross 
section. The critical and maximum Mach numbers are M,- and 
M,m respectively. A loss of about 5% occurs midway between 
M,- and M,m after which the loss increases considerably. Fig. 8 
illustrates the effect of Mach numbers on the pressure rise of a 
cascade of blades. 


Design Method. 


Certain generalised design curves have been evolved, and are 
helpful to the designer in fixing certain arbitrary figures. 


M, (37) 
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Example. 


Design an axial flow compressor given the following data :— 


Q = 50 lb. per sec. 

T, = 288°K. 

P, = 14-7 Ib. per sq. in. 
R.P.M. = 12,000. 

Q = 0:86. 

No. of Stages = 12. 

ATs = 35°C. 

Type of blading = 50% Reaction. 
Va = 600 ft. per sec. 


Tip diameter of compressor = 18 ins. 
The entry annulus*is determined from the equation of continuity. 


QO. P, 14-7 x 144 
A = 7h = ae = -_e_— 
Vino °° @! RxT, 2 ~  96x288 


= 0-0765 Ib. per cu. it. 


6 


0 O25 05 O75 
My 


Fig. 8—Effect of Mach. Number MI on the Pressure Rise in a Cascade 
of Compressor Blades. 
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50 
—_—_—_——. 1-090 sq. ft. 
600 x 0-0765 “ae tt 
1-090 x 144 
= 2 _ @®) from whi 4) 2 
Also A a/4 (D ) from which d@ 07a 18 
= 11-180. 
18-000 — 11-180 6-820 . 
Height of blade at entry = ee er es = 3-410 ins. 
ATs 


$2 0-923 (U/100) (V,/100) 
where U is velocity. 
12,000 . 7 . 14-590 


At the mean radius U = ea eee 763 ft./sec. 


_ 25-88 
0-923 . 7-63 . 6 . 0-86 
For 50% reaction , = $x Va/U (tan oc, + tan o,) (38) 


From equation (20) (tan oc,—tan oy) = 


(tan oc, -tan &,) = 0712. 


600 ; 
05 = 4. 763 (tan oc, +tan oz). 


Whence (tan oc, +tan oc.) = —— = 1-272. 


tan oc, — tan oO, = 0:712. 
tan oc; + tan oO, 1-272. 


ll 


Subtracting -2tan «<, = —0-560 
Adding 2 tan oc, = 1-984 


1-984 
2 


— 0-560 
tan Lo = a 0-280 Pa Le = 15° 39’ 


-2 
Since there is 50% reaction oc o, = 15° 39’ 
©, = Oo, = 44° 47° 
These angles occur at the mean radius r,, and it is now necess 
to obtain variations of these angles over the total blade height. 
Assuming that axial velocity and work input per lb. are constant 
with radius 
Va = constant (K,) (39) 
D 


tan oc, = 0992 .. oc, = 44° 47’ 


I 


I 
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and work done per Ib. of air per sec., W = U/g (Vy — Vwa) ft. lb. (40) 
where Vwa and Vw3 are the whirl components of the absolute 
velocities of the air at the inlet to and exit from the rotor. 


From equation (40) U, (Vw3— Vwo) = U7. Va (tan ocy~tan oc,) 


= constant (K,) (41) 
Where Ur = peripheral velocity at any radius r = U xr (42) 
In order to simplify the calculation 7 = 1. 
r 
Now Va 7 tan o, + tan o, = tan oO, + tan Oo, (43) 
From equations (41) and (42) U.7.K, (tan o,-tan «,) = K,, 
- K, 
or (tan oc,-tan <,) = ——— 
( $ o) U.7.K, a4) 
. U.r 
From equation (43) tano, + tana, = K (45) 
1 
; U.r 
and tan o, + tana, = (46) 
K, 
constant 
From the theory of free vortex V,. = — 
constant 
Va tana, = —— 
y 
Kg 
hence tano, = >. 
Substituting this in equations (44) and (45) we get :— 
' _ K2 n K3 
an %3 ~~ U.r.K, r 
U.r K3 
andtano, = =" - — 
K, Y 
assuming 50% reaction at the design radius 7 = 1. 
tan oc, = tan og when 7 = 1. 
K, U 
+ K -~ K 
U.K, $ K, 3 
U- Ky 
From which tan o, = —__U (47) 


iy 
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K, 
‘ — Ur | hae U 
ae K, 2K,.7 i (48) 
K 
U = 2 
tan oc Ke uy 49 
2 > Ul7.K, OK, .+ (49) 
K, 
‘ U.r ie U 
pee eo 2K, .7 0) 
For constant 50% reaction at all radii :— 
t tan oc, = 3 Cee Ke 5 
= n = ey 
le ne K, | U.K.? to) 
U.r7 K, 
tan oc, = tan oc, = 4 (52) 


K, 7 U.K,.r7 


i] 


763-600 (0-712) 
325-953. 


We can now obtain the air angles 0¢,, oc, and oo, oc, over the 
entire annulus. The air angle distributions are shown graphically 
in Fig. 9. 

From equation (11) « * = 0-8 (29° 08’) = 23° 18’ and from 
Design 
Deflection Curves S/C = 1-1 (see Fig. 10) obtained from cascade tests. 


[023 +01 | VJT1. 6 
0-275 6 | 
Now @ B,-B, and B, = «,-8 
] By -—(oc.+8) = B,—- c,+0-275 6 
¢ = ,-f, = 2:-2° (obtained from cascade) 
By = 44° 47’-2° 12’ = 42° 35’ 
42° 35’ — oc, +0-275 6 
0:725 8 = 42° 35’—15° 39’ = 26° 56’ 
@ = 37° 45’ and B, = B,-6 = 4° 50’ 
Also £ = f,—-6/2 = 42° 35'—18° 52’ = 23° 43° 


K, U.Va (tan oc, -tan «,) 


fr = Ove 


Substituting in equation (12) 6 


I 


D 
fl 
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Fig. 9—Air Angle Distributions—Constant Reaction. 
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Assuming that there are 42 blades in a row, at 7m :-— 


2 x 7:295 
S = a9 = 1-090 
Ge Ey eS 0p 
1-1 1-1 
Aspect ratio = ah = oa. = 3-450 
c 0-990 


This process has now to be repeated at selected points along the 
blade length. The S/C ratio is derived from the corresponding 
air angles and the chordal length is determined thus a complete 
picture of the blade is obtainable. The following table gives the 
complete data for the stage in question. 


ce or 
43° 36° 
44° 47°] 15° 39° 1-090 


46° 38’ | 25° 1° 1*309 


Adopting a Gottingen 436 aerofoil section, Fig. 11 illustrates 
the velocity diagrams and blade layout for the mean radius position. 


Drag and Lift Co-officients. 


From equation (6) tan o, = } (tan 44° 47’ + tan 15° 39’) = 

} (09924 +0-2801) = 4. 1-2725 = 0:6362 .°. om = 32° 28’, 

From equation (10) C, = 2.1-1 (0-7123) 0-8437 = 2-2 .0-7123. 
0-8437 = 1-325. i 

Hence from equation (26) C,, = 0-018 (1-325)? = 0-0315 

At this radius S = 1:090 and & = 3-410 .-. from equation (27) 

1-090 . 
3-410 
The variation of the profile drag co-efficient C,, is comparatively 


narrow for a wide range of geometrical forms of cascade and a 
constant value of 0-018 has ‘been adopted. 


C, = C.p4+C,.+C,, = 0:018 +0-031+0-006 = 0-055 


a5 = C, + sf Em 
°*] 


C,, = 0-020. = 0-006 


From equation (28) ae 
os* oc, 


24 


% 


AIR DEFLECTION € DEGREES 
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AIR OUTLET ANGLE 2° 


Fig. 10—Design Deflection Curves. 
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_ 0-055 cos? 19° 36’ 09-0506 
eo cos? 44° 47" 
From equation (3) the theoretical pressure rise through the 
blade row is found by putting the loss w equal to zero, which gives :— 


APth = 34 @ Va? (tan? oc,—tan? o,) = 
4 @ Va? (sec? oc, —sec? 0). 


APth ; sec?a, APth : cos?cc, 
oo, OE LE TOC = l1- —,.~ 
4@ Va? sec? oc, sec?oc, s0V,? COS? Ly 
0-503 
= -~ — = 0-458. 
1 ~ "0-927 Se 
0-0596 
From equation (32) 7, = 1 - Odea. 1-0-1830 = 87%. 
Substituting static temperatures in equation (21) we get 
7S . aa ( 0-870 . 25-88 4 3-3 
= (140) = (1 + 
a ( To 2 ms (CS 


= (1-0782)*5 = 1-3015. 


Overall Performance. 


The overall static pressure ratio R = = = (yas 
P, T, 
s (3) ra (53) 
Ty 


We have assumed that T, = T, where T, = the ambient temper- 
ature, but in fact it will be less than T, due to the temperature 
equivalent of the air velocity at the mean radius. This velocity 
is V, (see Fig. 7). 

(Vasecoc,)? (600 sec 15° 39’)? 


Temperature equivalent of V, = ——_~> = ——— 
' i ° Qe. J.Cp 21,638 


17-9°C. 


T, = 288-179 = 270-1°K. 
0-87 . 25-88) °° 
R, = ( eae = (1-0834)85 = 1-336. 
Assuming Va, is constant throughout the compressor 
T,-T, = 12x25°88 = 310-56°K. 
a. Te 310-56 +270-10 = 580-66°K. 


26 COMPRESSOR DESIGN FOR AIRCRAFT ENGINES 


Fig. 11—Layout of Compressor Blading. 
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580-66 °°" (sans me 
) z-10) 


| 


P 
From equation (53) = rer 
1 


10-283. 


l 


Compressibility Effects. 
The relative velocity V, at the mean radius at the first stage is 
given by :— 
V, = Vasec ©, = 6001-4088 = 845-28 ft./sec. 
M 845 845 
~ af B22x14x96x27010  —«:1170 
This figure is quite satisfactory. 
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CENTRIFUGAL COMPRESSOR DESIGN 
FOR AIRCRAFT ENGINES. 


The rapid development of gas turbines during the war resulted 
in the majority of aero-engine manufacturers adopting the centri- 
fugal compressor rather than its axial counterpart. Previous 
design experience with superchargers, based on centrifugal principles 
no doubt influenced this trend. High delivery pressures were 
initially not required as total thrust was not of such a high order 
as it is to-day and weight was an important factor as was total 
length. All these factors resulted in the centrifugal compressor 
being favoured as a means to secure adequate compression. 

Three types of impeller vanes are possible, namely, the back 
swept vane, the straight radial vane and the forward swept vane. 
The three types are illustrated in Fig. 1. The back swept and 
straight radial types are in common use, but the forward swept 
type is not used, despite the fact that it gives the highest pressure- 
ratio. Its chief drawback lies in the excessive diffuser losses, due 
to the considerable absolute air velocity leaving the impeller. 

In comparing the straight radial with the backward curved 
type it can be seen that the straight radial vane is a particular case 
of the curved backward vane in which the blade entry and outlet 
angles are equal. Therefore any calculations for one type are 
applicable to the other, provided allowances are made for the 
differences in vane angles. 


Nomenclature. 
Va, = Absolute velocity of entry air in feet per second. 
V, = Velocity of entry air relative to impeller in ft./sec. 
Vwo = Whirl velocity in feet per second. 
U._ = Tangential velocity of impeller tip in feet per second. 
V. = Absolute velocity of outlet air in feet per second. 
U, = Velocity of outlet air relative to impeller in ft./sec. 
U, = Tangential velocity of blade edges at entry in ft./sec. 
w = Angular velocity of impeller in radians per second. 
Vio = Absolute velocity of entry air with prewhirl in ft./sec. 
Vig = Radial velocity in feet per second. 
o = Slip factor. 
y, = Radius of entry annulus in feet. 
v, = Radius of impeller in feet. 
p = Power input factor. 
Q = Mass flow of air. 
T,: = Inlet total-head temperature. 


Inlet total-head pressure. 


= 
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Tot = Outlet total-head temperature. 
= Outlet total head pressure. 


at 
y = Ratio of specific heats 2 
ny; = Isentropic efficiency of compressor. 
oc, = Vane inlet angle. 
C, = Specific heat at constant pressure. 
C, = Specific heat at constant volume. 


Compressors can be single or double-sided as shown in Fig. 2. 
(a) and (8) or should higher compression ratios be required they 


. SSSoSosony 


Ve 


Fig. 2—Centrifugal Compressors. 


"cao a ceptor ad ncaa 
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can be of the compounded type. Basically their fundamental 
theory is the same. The majority of modern centrifugal com- 
pressors are of the straight radial type. 


Air is sucked into the impeller eye at an absolute velocity Va, 
and a relative velocity V, is whirled at a peripheral speed U and 
flung out by centrifugal force at an absolute velocity V, and relative 
velocity U,. The centripetal acceleration is obtained by a pressure 
head, therefore the static pressure of the air increases from eye to 
tip of impeller. The kinetic energy of the air is now converted 
into pressure energy by means of a diffuser by which the final 
velocity of the air is reduced to somewhere in the region of the 


- velocity with which the air originally entered the impeller eye. 


Work Done During Compression. 


Fig. 3 illustrates a typical centrifugal compressor with the 
necessary nomenclature employed. It can be seen that the 
incoming air has an initial velocity of V., and that it possesses no 
angular momentum. In order that the air can pass smoothly into 
the vanes, the axial portion between the eye root and eye tip must 
be curved. This inlet angle o, will be given by the direction of 
the relative velocity V,. 


Under the action of centrifugal force the air leaves the tip of 
the impeller with an absolute velocity V, and a tangential or whirl 
velocity Vw, together with a small radial component V,,. Under 
the ideal conditions the velocity V, would be such that the whirl 
velocity Vw = U, the tangential speed of the impeller tip. How- 
ever, due to the reluctance of the air trapped between the vanes 
to move round, with the impeller, the static pressure on the leading 
face of the vane exceeds that of the trailing edge so that the air is 
prevented from acquiring a whirl velocity equal to the impeller 
speed. This phenomenon is called “‘slip” and its amount is largely 
dependent upon the number of vanes in the impeller. In practice 
this figure is smaller when the number of vanes are increased and 
it is necessary in the course of design to assume a value for ao where 


Vwe 
is the ratio —— 
. U 


Considering a flow of 1 1b. of air per second through the impeller 
the theoretical torque is equal to the rate of change of angular 
momentum and is given by :— 


Ve. 7 
Theoretical torque = a Ib. ft. (1) 


If the angular velocity is radians per second is w the work done 
on 1 Ib. of air per second will be :— 


33 
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Vws X%o XW Vwe x U 
a OC ft. Ib. per sec. (2) 
§ § 
When introducing the slip factor, theoretical work done 
o Ue 


ft. lb. per second. 


During the passage of the air through the vanes, certain losses 
occur such as friction, etc., which result in the actual work input 
being greater than the theoretical value given in (2). To com- 
pensate for this a “power input factor’’ is introduced so that the 
actual work done on the air becomes 

ees ft. lb. per second (3) 
§ 

If Q. lb. of air per second are passing through the impeller then 

the total work done per second by the impeller is 


U2 
Sp SE & ie. fiee, 
§ 
and the horse-power equivalent of the work done on the air is 
4 x U? 
fj. ( 
530 g 


Pressure Rise. 


The work done on the air during compression can be considered 
in terms of temperature changes across the whole compressor. 
Ideally the type of compression should be isentropic (7.e., reversible 
adiabatic) but due to the friction of the air passing through the 
vane channels, a perfect isentropic compression cannot be obtained 
hence the introduction of the power input factor. 


The difference between the power input factor and the slip 
factor should be clearly understood. Introducing the power 
input factor represents the increase in total work input due to 
overcoming frictional losses which results in the raising of the 
outlet temperature. It then follows that the maximum cycle 
temperature can be reached without burning so much fuel, hence 
these losses are not entirely wasted. Thus we can still adopt the 
isentropic type of compression for calculation purposes, bearing 
in mind that the power input factor should be kept as low as possible. 
Present day designs require that this factor should be between 
1-030 to 1-040. 


_ © 
The slip factor is the ratio > and it can be seen from the 


velocity triangles in Fig. 3 and its introduction into formulae 
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determining total work decreases this figure. It will be clearly 
seen that the nearer Vw, reaches the value U the greater the total 
work done. Increasing the number of vanes in the impeller would 
appear to increase o, but in doing so the effective flow area of the 
impeller eye is decreased. Greater frictional losses thus occur, 
since for the same mass flow the inlet velocity must be increased. 
A compromise has been found by experimental means whereby a 
suitable number of vanes will give a slip factor of 0-9, 7.¢., about 
25 vanes. 


Now the temperature equivalent of the work done on the air 
given by equation (3) 

xa x U? 2 

- 2xex— (5) 
ExJx& 

where Cp is the mean specific heat over the temperature range 


Tot — Ty, which is the total-head temperature rise across the 
impeller. It therefore follows that :— 


pxaxU? 

Tor —T — 6 

( at 11) g x J xC ‘bp ( ) 

Also work done = Q J Cd (Ta —Tyt) ft. Ib./sec. (7) 
Horse-power equivalent = DeLee (Tot — Tt) (8) 


Adiabatic Efficiency. 


If a value be assumed for the total-head isentropic efficiency, 
nc, then we know how much of the work is actually used to compress 
i -air. Thus the total-head pressure ratio may be found as 
follows :— ; 


as aa 
Poo ae) ry, tT 
Pio|}|©60CiaTQtC [ sais aaa 
me 
= [iene » ie LE i (9) 
gxJxCpxTi 


From equation (9) it will be seen that the factors influencing 
the pressure ratio are the power input factor #, the slip factor a, 
the impeller tip speed U and the total-head inlet temperature T,?. 
A well-designed centrifugal compressor when running under near 
perfect conditions can achieve an efficiency of 75% -80%. Very 
little energy is actually lost during compression but quite a con- 
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siderable portion of the energy goes into heating the air to a temper- 
ature higher than that which would be obtained from a perfect 
isentropic compression. 


The ideal temperature rise corresponds to a perfect isentropic 
compression, 1.¢., T,¢1-T,#, and if the actual total head temperature 
tise is T,¢—T,¢ then it follows that :— 


P-Ty 
nc = Tye Te (10) 
Tit -Ty 
yr 
Paty ¥ 
oy L ipa) = ] 7 
Be ra Tits (11) 


Sometimes it is convenient to express the isentropic efficiency 
in terms of pressure-ratio, which can be done in the following 
manner :— 


Work required for perfect isentropic compression is given by 
the expression :— 


QxJxCp (T,t-T,2) ft. Ibs. per second. 


and the actual work required to effect the compression is given by 
the expression :— 


QxJxCp (T,¢-T,Z) ft. lbs. per second. 


yru 
Pas: 

But T, = Ty (Si) 
a4 

& P as n 
and Tj = Ty, ea 


Substituting these values-into expressions for work done :— 
yo1 


P,t 
Isentropic work = QxJxCpxT;t [ re, -1 ] 
1 


Actual work 


QxJxCpxTy EG -1 | 
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Since the isentropic efficiency is the ratio :— 


isentropic work 
actual work 


then nc = (24) * = 
(12) 


; 2. oe , 
Denoting the pressure-ratio ¢ by 7 the expression (12) becomes 


1 


y-1 
yo 
sp tes ee (13) 
D-1 
i | 


Theoretically, it is very difficult to obtain the value of 7 and 
therefore recourse should be made to the data obtained from tests 
carried out on previous compressors. 


: We will now consider the design of a double sided cones 
compressor having been given the following data :— 


Overall diameter of impeller 21 ins. 
Eye tip diameter 12 ins. 
Eye root diameter 6-5 ins. 
Total mass flow Q 40 Ib. per second 
Power input factor p 1-040 
Slip factor ¢ 0-90 
Revs. per minute 15,000 
Inlet total temperature T, 288°K . 
Inlet total pressure P, 14-7 lbs. per sq. in. 
Total head isentropic ‘efficiency 76% 
Density of inlet air 9 _ 0-0765 Ib./cu. ft. 
15000 x 21 xz 
U 60 x 12 1371 ft. per sec. 
1:04 x 0-9 x 1371? : 
From (6) T,f-T,t = 322x1400x004 7 162-6°C. 
Toty-Tit = 0°76x1626 = 123-6°C. 


40 x 1400 x 0-24 x 162-6 


8) H.P. ired 
From (8) require 550 


= 3973 H.P. 
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The inlet velocity which is axial is not known and in order to 
determine the inlet angle of the vanes it is necessary to satisfy 
the equation of continuity Q = 9, x Ax Va,, and it is assumed that 
the mass flow is equally divided between both sides of the impeller. 


a (12? - 6-52) 
A = eis) = 0:555 
rs ies 0-555 sq. ft. 
The equation of continuity Q = AxVa,x@ as 9 = 0-0765. 
20 = 0:0765 x 0-0555 x V., whence V., a 
0-0765 x 0-555 
= 469-4 ft. per sec. 
The inlet total head temperature T,, = T, + — 
2gxJxCp 
= per. 
64-4 x 1400 x 0-24 
= 288°K + 10-2°C = 298-2°K 
y 
op, Ti \yva 298-2 \** 
Now eam 2 — ce. H = A es 
P, T, ence .P,; 14-7 x 988 
= 14:7x1:125 = 16:5 lb. per sq. in. 
Pt 123-6 7° 
From (9 i E + ] = 3 
(9) Pt 298-2. sia 
| ee 7x1 
The peripheral speed at the eye tip radius = _ 
= 786 ft./sec. 


The peripheral speed at the eye root radius 


~ PERORINO _ Jo9 i. per sec. 


12-60 
469 
co, at root = tan-} 195 7 47° 49’ 
; 469 
oc, at tip = tan} 786 30° 49’ 


The Diffuser. 


The objects of the diffuser are :— 

(1) To provide a smooth canalization of the flow of air from 
the impeller and to reduce to a minimum any turbulence which 
would otherwise be set up in the volute casing or ducts. 
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(2) To convert the kinetic energy to pressure energy. It must 
be stressed that it is much more difficult to obtain an efficient 
declaration of flow, than it is to obtain an efficient acceleration. 


With the compressor running under its optimum conditions 
the vanes of the diffuser can be designed to pick up the air smoothly 
but there is a natural tendency for the air to break away from the 
surface in a diffusing system thus forming areas of turbulence. 
Consequently some of the kinetic energy is transferred into thermal 
energy instead of useful pressure energy. Undesirable temperature 
rises are also caused by the frictional effect of the high speed fluid 
flow through the diffuser channels. Experimental work has 
shown that the maximum permissible divergence of the diffuser 
channel is about 11°. Any increase in this figure results in higher 
losses. 

We will now consider what happens to the air on leaving the 
impeller if the diffuser were not there, Fig. 4. 


PATH OF AIR LEAVING IMPELLER 
TIP 15 EQUIVALENT TO A 
LOGARITHMIC CURVE 


Fig. 4—Air Flow after Impeller Tip. 


On leaving the impeller the air has an absolute velocity V., a 
tangential velocity Vw, and a relatively small velocity V7,. It 
follows that upon leaving the impeller its angular momentum 
remains constant, as no other forces act upon it which would alter 
its momentum. The angular momentum at any point is given by 
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x Vw,xr : 
DANees and as the mass and momentum is constant the tan- 


gential component Vw, of the absolute velocity must be inversely 
proportional to the radial distance. Therefore the radial velocity 
Vr, will also decrease from the impeller tip to the diffuser vane, in 
accordance with the continuity equation. Since both components 
are inversely proportional to the radius, the absolute velocity is 
also proportional to the radius. Furthermore the angle which 
the direction of the absolute velocity V, makes with the tangential 
velocity Vw, at any point is constant. Thus we have the con- 
ditions for a logarithmic spiral which is the theoretical path of the 
air leaving the, impeller. 

Both Vr, .and Vw, decrease after leaving the impeller tip then 
the resultant velocity V, will also decrease. This occurs in the 
vaneless space between the impeller tip and the diffuser and some 
diffusion will take place. For reasons which will be explained 
later the number of diffuser passages is less than the number of 
impeller blades. 


We will now consider the design of a diffuser to be used with 
the impeller designed earlier. Certain dimensions have to be 
assumed and are as follows :— 


Number of diffuser vanes, 13 

Width of casing, . 1-750 ins. 
Radial width of vaneless space, 1-500 ins. 
Mean radius of diffuser throat, 12-75 ins. 


It will also be assumed that any losses are equally divided 
between the impeller and the diffuser. From the above data, 
radius of vane leading edges = 12 ins. 


Vwr xv = constant. .*. whirl velocity of air at 12 ins. radius 
S. ge STRAWS: 1996-46 per sec. 
12 
1079? 1,164,241 e 
T . = = > , pe 3°8° 
emperature equivalent Dex] xCp "a ESL —21,638-4 rs) 


Now the radial velocity can be obtained from the continuity 
equation and as the density depends on the radial velocity Vr, a 
trial and error method must be used. The density changes only 
slightly for considerable changes in the radial velocity V7,, hence 
at least 5 figure logarithims should be used. Only the final trial 
will be given. 

Assume Vr, to be 290 ft. per sec. 

290? 


Temperature equivalent “= dex ]xcp = 3-8°C. 
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Temperature equivalent of resultant velocity = 53:8+3-°8 
57-6°C. 
460-8°K. 


Tt —T zt = 162-6°C. .. T,¢ = 298-2 + 1626 
But T,¢ is the same for any radius after the impeller tip. 
T, at 12” radius = 460-8-57:6 = 403-2°K. 

To find the static pressure P, and the density @, we require 
to know the static temperature, on the assumption of isentropic 
compression up to the 12 in. radius. 

“. Ty — Tot = 39°C. 
As 50% of the losses occur before the diffuser T,’ at the 12 in., 


radius = 403-2 — 1955 = 383-7°K. 
P T1 va 383-7 °° 
2 2 . . 
—2 _ (2 & 2in, radius = 165 ( ) 
Pt Geo P, at 12 in, radius 998-2 
= 1652-417 = 39-882 lb. per sq. in. 
P, = oe = 0-1483 Ib. per cu. ft. 
Area of cross section of flow at 12in. radius = rae 
= 0-916 sq. ft. 


Thus having Q, @g2 and A we can check on the radial velocity 
Vr2 by means of the continuity equation :— 


Vr2 = oe 294 ft. per sec. 
0-1483 x 0-916 Vr 
The angle of the diffuser vane leading edge is tan ( a ) 
2 
= tan“ (=) = 15°-15’ 
1079 


The next thing to consider is the throat width of the diffuser 
channels and this can be obtained by a similar process. The mean 
radius of the diffuser throat is 12-75 ins., and as the losses between 
a radius of 12 ins. and 12-75 ins. practically remain the same, we 
shall still consider 19-5°C. as this loss before the diffuser throat. 
Proceeding as before :— 


Whill'velocityet 1975 anecad, @: oe 


12°75 
= 1016 ft./sec. 
Temperature equivalent = 2 SOO ads 47-7°C. 
2gxJxCp 


Try a radial velocity of 275 ft. per sec. 
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. 2752 
Temperature equivalent = ~~~ = 3-5°C. 

7 p quivalen Dex J xp 3-5°C 
Temperature equivalent of resultant velocity =47-7 + 3-5 =51-2°C. 
T2 at 12°75 in. rad. = 460°8-—51-:2 = 409-6°K. 

T2! at 12°75 in. rad. = 409-6-19-5 = 390-1°K. 


3°5 


390-1 
Ph = 165 ( ) = 165x2560 = 423 Ib. . in, 
2 298-2 16-5 x 2-560 = 42-3 lb. per sq. in 
144 x 42-3 
Cc =~ = 0-155, lb. . ft. 
2 96 x4096 0-155, lb. per cu. ft 


Neglecting the thickness of the vanes the area of flow in a radial 
2m X 12-75 x 1-75 


directi: a ) ae =) S0-978 sa. fet. 
mn 144 “4 
Vr, at 12°75 radius _ 40 = 266 ft. per sec 
2 0-155 x0973 * Per See. 
Angle of direction of flow = tan™ 016 7 14° ~— 40 ft. 


The total throat area of the diffuser passages equals 0-973 sin. 
14° — 40’ = 0-246 sq. ft. 
The width of each throat must therefore be ciacas os 
xT 


= 1-55 ins. 

In order to obtain the correct curvature of the vanes it is 
necessary to determine the logarithmic spiral from the data already 
obtained. It would be possible to obtain successive values for V7 
and Vw,, from which we would obtain-a series of values for ¢ but 
this would be a lengthy process, hence we must resort to the 
logarithmic spiral. 

The equation to a logarithmic spiral is 7 = ae@ where 7 = ray 
or radius length ; 0 is the coefficient of @ =cotocwhere oc =90° —¢ ; 
§ = ray angle and a = initial radius. Obviously knowing ¢ at 
12 in. radius this is a convenient point at which to commence the 
logarithmic spiral. Hence oc = 90°-¢ = 90°-15°. 

= 75° and cot 75°. 
= 0-26794. 


0-807 040 
The equation now reads r = 12e 
. 12 0-26794a 0.004678 
rn. 


In log form log 7 log 12 + 0-00467a log e. 
1-0791812 + (0-00467 x 0-4343a). 


1-0791812 + 0-0020282. 


+ 
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1-16028 | 1-17044 


35° 


04056 | 0-05070 | 0-06084 | 0-07098 | 0-08112 | 0-09126 
13-805 | 14-131 | 14-464 | 14-806 


20° 


15° 
03042 | 0 


TABLE OF VALUES. 
0-01014 | 0-02028 | 0 


08932 1.09046 | 1-10960 1-11974 | 1-12988 | 1-14002 | 1-15016 
| 12-574 | 12-871 


12-283 


1-0791812 | 1- 


The above values can now be plotted on to the diffuser layout 
Fig. 5 where by using the appropriate radius to connect the plotted 
points we obtain the correct diffuser vane curvature. It will be 
seen from Fig. 5 that in order to obtain optimum conditions the 
direction of air flow after the throat should be in the direction as 
previously calculated, 7.e., 14° 37’, this being the angle of flow 
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at mean radius of throat area. However some amendment may 
be necessary in order to accommodate the diffuser into the general 
design of the engine but should be avoided as much as possible. 


Compressibility Effects. 


When the velocity of a compressible fluid, relative to the surface 
over which it is moving, reaches the speed of sound in the fluid, 
pressure losses occur. 

When designing a centrifugal compressor the size of the unit 
is an important feature and the largest possible mass flow with the 
smallest compressor is the optimum aim. Air speeds are high, 
thus it is of great importance that the Mach numbers at certain 
points in the compressor do not exceed the value beyond which 
losses occur dué to the formation of shock waves. Normally, 
the Mach numbers should not be allowed to exceed 0:8. 

We will now consider the Mach numbers at various stages in 
the compressor, 7.¢., (a) Impeller Intake; (6) Impeller Passages ; 
(c) Diffuser. 


(a) Impeller Intake. 


Due to the curvature of the impeller blades the air at entry is 
deflected through a certain angle before it passes into the radial 


~ channels. There always exists a tendency for the air to break 


away from the trailing edge of the impeller eye, hence here is a 
point where the Mach number is extremely important and is shown 
in Fig. 6 (a). 

Fig. 6 (b) illustrates the velocity triangle for the impeller eye 
and is shown in full line. It has been assumed that the inlet 
velocity Va, is constant all over the surface of the impeller, therefore 
the velocity of the entry air relative to the impeller designated 
by V, will reach a maximum at the eye tip of the impeller. The 
inlet Mach number will be given by :— 

M = gota ies (14) 
Vv gyRT, 
where T, is the static temperature at the inlet. 

We will now consider the Mach number at inlet for the com- 
pressor being designed. 

From previous calculations :— 
Inlet velocity Va, 
Eye tip speed 
Relative velocity at tip V, 
Velocity of sound 


469 ft. per sec. 
786 ft. per sec. 
- 469? +786? =915 ft. per sec. 
V32-2 x 1-4 x96 x 288 — 
1116 ft. per sec. 


t 


ll 


to 
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Maximum Mach number at inlet = Ele = 0-82. 


1116 
This figure would be satisfactory under sea level conditions, 
but when operating at high altitudes with ambient temperature 
around 215°K. the Mach number would be much higher. 


Thus we will consider the conditions experienced at an altitude 
of 40,000 ft. The effect of drop in atmospheric temperature will 
not be so great as may be imagined due to the ram effect at the 
intake when the aircraft has a forward speed. Assuming a normal 
speed of 375 m.p.h. the following data is obtained :— 


550? 
2gxJxCp 
=: 13-98°C. 
Inlet total head temperature = 215°K +13-98°C = 228-98°K. 
Temperature equivalent of axial inlet velocity = 10-2°C. 
Inlet static temperature at 40,000 ft. = 228-98 — 10-20. 
= 218:°78°K. 


Mach number at inlet = 0:82 J an = 0-940. 


218-78 

This then exceeds the permissible Mach number and it is, 
therefore, necessary to reduce the relative velocity V, ‘It is 
possible to do this by means of introducing ‘‘prewhirl’ at the 
intake. Curved vanes attached to the compressor casing deflect 
the incoming air and we obtain a velocity triangle as shown in 
Fig. 6 (0). 

Unfortunately this method of reducing the relative inlet velocity 
V, reduces the work capacity of the compressor, by reason of the 
inlet air possessing a whirl velocity V\,. The rate of change of 
angular momentum per Ib. of air per second now becomes :— 


Temperature equivalent of forward speed = 


PREWAIAL ANGLE 
Ue 


a . 
¥ A 
INLET GUIDE VANE 


(a) (b) 
Fig. 6—-Effect of Prewhirl. 
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& 

If it be assumed that V.,, remains constant over the eye of the 
impeller, then the angular momentum must increase from root to 
tip. It follows that the amount of work done per Jb. of air varies 
according to the radius at which the air enters the impeller eye. 
However, the Mach number is only high at the eye tip and it is 
preferable to vary the prewhirl from root to tip by means of twisting 
the inlet guide vanes so that the prewhirl value varies from zero 
to a maximum at the eye tip. 


The angle of prewhirl is generally made 30° and a new value 
for the axial velocity must now be found by the usual trial and 
error process. 


Only the final result will be shown. 
Try Va, = 490 ft./sec. 
490 
cos 30° 
Temperature equivalent of Vz, = 14-75°C. 


Ib. ft. (15) 


Vag ™ = 565 ft./sec.. 


T, = 2982-14-75 = 283-45°K. 
FP ie 13-818 Ib i 
1 ~~ (1-052)3°5 a ° per sq. in. 
13-818 x 144 


0-0731 Ib. per cu. ft. 


el = 96%28345 


20 


7 (12? 6-52) 


Whirl velocity at inlet Vw, = 490 tan 30° = 282 ft. per sec. 
Maximum relative velocity = +~/490%+ (786 —282)? = 1/494,116 


Check on Vay = = 493 ft. per sec. 


= 703 ft. per sec. 
Maximum inlet Mac number when T,, = 298-2°K. 
703 
=... =_ 0-634. ? 


4/322 x 1-4 x96 x 283-45 


This is a satisfactory figure under altitude conditions. 


Prewhirl has an effect on the pressure ratio and assuming that 
the prewhirl remains constant over the whole of the impeller eye :— 


78 
as = 605-5 ft. per sec. 


Speed of eye at mean radius Ue = 5 
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Actual temperature rise _P (0 U?—Vous x Ue) 


gExj]xCp 
1-04 (0-9 x 13712 — 282 x 605-5) 
eee eee eee ss «146 -2°C, 
32-2 x 1400 x 0-24 6 
Tot, —Tit = 0-76 x 146-2 = 111-11°C. 
Pat 111-11 7°° 9-5 
= = ———— = (1-372 
Pit [ + "998-2 (1-372) 
= 3-02. 


Mach Number in the Diffuser. 


The maximum: Mach number in the diffuser obviously occurs 
at the entry, that is, the impeller tip. 


Vwe at impeller tip = 0-9x1371 = 1233 ft. per sec. 
Try Vi, = 390 ft. per sec. 
‘Resultant velocity V, = ~/1233?+390? = 1293 ft. per sec. 
Temperature equivalent = ie pep = 77-28°C. 
Tet = 460-8°K. .:. T, = 4608-77-28 383-52°K. 
Again assuming a loss of 19:5°C. up to this point. 

T,! = 3835-195 = 364°K. 


I 


Vv 


P T,1\ 7-2 364 \*> 
Now =e = (2) P, = 165 ( ) 
Pit Tyt ° > \ 298-2 
= 33:°148 Ib. per sq. in. 
33-148 x 144 
2 = ———— + = 0-129 lb. . ft. 
@ 96 x 383-52 9 Ib. per cu. ft 


Area of flow in radial direction at 10-5 rad. = 2x 105 x 175 


144 
= 0-801 sq. ft. 
Checking value of V,, = Se wT = 387 ft. per sec. 
-129 x0- 


Thus the resultant velocity of 1293 ft. per sec. is satisfactory 
and therefore the Mach number at the impeller tip is :— 


1293 


ey = 1 ‘003. 
4/32-2 x 1-4 x 96 x 383-5 
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We will now consider the Mach number at the diffuser leading 
tip edges. Using data already obtained the resultant velocity 


= 4/2902-— 10792 = 1117 ft. per sec. 
The static temperature T2 at this radius was 403-2°K. 
the Mach number is 
1117 
V/32-2 x 1-4 x 96 x 403-2 

From previous experimental work it has been found that Mach- 
numbers greater than 1-00 can be used at the impeller tip but that 
high Mach numbers present at the leading edges of the diffuser 
are undesirable. Thus can now be easily seen the reason for the 
vaneless space. Circumferential pressure variations are set up. 
when Mach numbers are high at the diffuser tip. These pressure 
variations cause vane vibrations which should they equal the 
natural frequency of the impeller then fatigue will occur. In 


order to reduce resonance the number of vanes in the impeller 
should not be a multiple of the vanes in the diffuser. 


= 0-84. 
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List of A.E.S.D. Data Sheets 


Safe Load on Machine-Cut Spur Gears. 

Deflection of Shafts and Beams. \ Connected. 
Deflection of Shafts and Beams (Instruction Sheet). 

Steam Radiation Heating Chart. ; 
Horse-Power of Leather Belts, etc. 

Automobile Brakes (Axle Brakes). \ . 
Automobile Brakes (Transmission Brakes). Connected 
Capacities of Bucket Elevators. 

Valley Angle Chart for Hoppers and Chutes. . 
Shafts up to 5} inch diameter, subjected to Twisting and Combined 
Bending and Twisting. . 
Shafts, 53 to 26 inch diameter, subjected to Twisting and Combined 

Bending and Twisting. 
Ship Derrick Booms. 
Spiral Springs (Diameter of Round or Square Wire). 
Spiral Springs (Compression). 
Automobile Clutches (Cone Clutches). 
” » (Plate Clutches). 
Coil Friction for Belts, etc. 
Internal Expanding Brakes. Self-Balancing Brake 
hoes (Force Diagram). 
Internal Expanding Brakes. Angular Proportions Connected 
for Self-Balancing. 
Referred Mean Pressure Cut-Off, etc. 
Particulars for.Balata Belt Drives. 
#” Square Duralumin Tubes as Struts. 
|’ 


3” Square Steel Tubes as Struts (30 ton yield). 


$ » » » (30 ton yield). 
1” 0 ” » (30 ton yield). 
a » ” ” (40 ton yield). 
’ » ” ” (40 ton yield). 
1” » (40 ton yield). 


Moments of Inertia of Built-up Sections (Tables). 
Moments of Inertia of Built-up Sections (Instructions Connected. 
and Examples). . 
Capacity and Speed Chart for Troughed Band Conveyors. 
Screw Propeller Design (Sheet 1, Diameter Chart). 
” ” ” (Sheet 2, Pitch Chart). Connected. 
” ” » (Sheet 3, Notes & Examples). 
Open Coil Conical Springs. 
Close Coil Conical Springs. 
Trajectory Described by Belt Conveyors (Revised, 1949). 
Metric Equivalents. , 
Useful Conversion Factors. 
Torsion of Non-Circular Shafts. 
Railway Vehicles on Curves. 
Coned Plate Development. 
Solution of Triangles (Sheet 1, Right Angles). 
Solution of Triangles (Sheet 2, Oblique Angles). 
Relation between Length, Linear Movement and Angular Movement 
of Lever (Diagram and Notes). 
» ” a? : ” » ” (Chart). 
Helix Angle and Efficiency of Screws and Worms. 
Approximate Radius of Gyration of Various Sections. 


92. 


Helical Spring Graphs (Round Wire). 
it % rr (Round Wire). Connected. 
me 53 5 (Square Wire). 
Relative Value of Welds to Rivets. ; . 
Graphs for Strength of Rectangular Flat Plates of Uniform Thickness. 
Fy Deflection = re ” ” ” 
Deflection of Leaf Spring. 
Strength of Leaf Spring. 
Chart Showing Relationship of Various Hardness Tests. 
Shaft Horse-Power and Proportions of Worm Gear. 
Ring with Uniform Internal Load (Tangential Strain) Connected 
sak’ 7 7 (Tangential Stress) . 
Hub Pressed on to Steel Shaft. (Maximum Tangential Stress at Bore 
of Hub). 
Hub Pressed on to Steel Shaft. (Radial Gripping Pressure between 
Hub and Shait). 
Rotating Disc (Steel) Tangential Strain. 
95 z i Stress. 
Ring with Uniform External Load, Tangential eens 
PY 5 z5 i 3 tress. 
Viscosity Temperature Chart for Converting Commercial 
to Absolute Viscosities. 
Journal Friction on Bearings. Connected, 
Ring Oil Bearings. 
Shearing and Bearing Values for High Tensile Structural Steel Shop 
Rivets, in accordance with B.S.S. ga i i 
Velocity of Flow in Pipes for a Given Delivery. 
Delivery of Water in Bipes for a Given Head. \ Connected: 
(See No. 105). 
Involute Toothed Gearing Chart. 
Variation of Suction Lift and Temperature for Centrifugal Pumps. 
Curve Relating Natural Frequency and Deflection. 
Vibration Transmissibility Curved or Elastic Suspension. Cc ted 
Instructions and Examples in the Use of Data Sheets, OREEGSS 
Nos. 89 and 90. 
Pressure on Sides of Bunker. 


Connected. 


Connected. 


93-4-5-6-7. Rolled Steel Sections. 
98-9-100. Boiler Safety Valves. 


102. 
103. 
104. 
105. 
106. 


Pressure Required for Blanking and Piercing. 

Punch and Die Clearances for Blanking and Piercing. 

Nomograph for Valley Angles of Hoppers and Chutes, 

Permissible Working Stresses in Mild Steel Struts with B.S. 449, 1948. 

Compound Cylinder (Similar Material) Radial Pressure of Common 
Diameter (D1). 


(Data Sheets are 3d to Members, 6d to others, post free). 


Orders for Pamphlets and Data Sheets to be sent to the Editor, 
The Dvraughtsman, cheques and orders being crossed ‘‘A.E,S.D.” 


